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Abstract. From June to July 2009 more than thirty differ- 
ent in-situ and remote sensing instruments from all over the 
world participated in the Cabauw Intercomparison campaign 
for Nitrogen Dioxide measuring Instruments (CINDI). The 
campaign took place at KNMI’s Cabauw Experimental Site 
for Atmospheric Research (CESAR) in the Netherlands. Its 
main objectives were to determine the accuracy of state-of- 
the-art ground-based measurement techniques for the detec- 
tion of atmospheric nitrogen dioxide (both in-situ and re- 
mote sensing), and to investigate their usability in satellite 
data validation. The expected outcomes are recommenda- 
tions regarding the operation and calibration of such instru- 
ments, retrieval settings, and observation strategies for the 
use in ground-based networks for air quality monitoring and 
satellite data validation. Twenty-four optical spectrometers 
participated in the campaign, of which twenty-one had the 
capability to scan different elevation angles consecutively, 
the so-called Multi-axis DOAS systems, thereby collecting 
vertical profile information, in particular for nitrogen diox- 
ide and aerosol. Various in-situ samplers and lidar instru- 
ments simultaneously characterized the variability of atmo- 
spheric trace gases and the physical properties of aerosol 
particles. A large data set of continuous measurements of 
these atmospheric constituents has been collected under var- 
ious meteorological conditions and air pollution levels. To- 
gether with the permanent measurement capability at the CE- 
SAR site characterizing the meteorological state of the atmo- 
sphere, the CINDI campaign provided a comprehensive ob- 
servational data set of atmospheric constituents in a highly 
polluted region of the world during summertime. First de- 
tailed comparisons performed with the CINDI data show 
that slant column measurements of NO 2 , O 4 and HCHO 
with MAX-DOAS agree within 5 to 15%, vertical profiles 
of NO 2 derived from several independent instruments agree 
within 25 % of one another, and MAX-DOAS aerosol optical 
thickness agrees within 20-30% with AERONET data. For 
the in-situ NO 2 instrument using a molybdenum converter, 
a bias was found as large as 5 ppbv during day time, when 
compared to the other in-situ instruments using photolytic 
converters. 


1 Introduction 

The Cabauw Intercomparison campaign for Nitrogen Diox- 
ide Measuring Instruments (CINDI) was conducted in June- 
July 2009 with the purpose of intercomparing and intercal- 
ibrating ground-based instruments measuring nitrogen diox- 
ide (NO 2 ), of which the observations can be used for the val- 
idation of tropospheric NO 2 retrievals from satellite observa- 
tions. The campaign was based in Cabauw, the Netherlands, 
a place with mainly agriculture in the direct proximity but 
at <40 km distance from the 4 largest cities of the Nether- 
lands, the closest of which is Utrecht (18 km). The Cabauw 


Experimental Site for Atmospheric Research (CESAR) is 
managed by the Royal Netherlands Meteorological Institute 
(KNMI). CINDI was building on previous and more ex- 
ploratory campaigns focusing on air pollution measurement 
techniques and satellite data validation (e.g. DANDELIONS 
2005 and 2006 in Brinksma et al., 2008; Hains et al., 2010, 
and FORMAT in 2002 and 2003, see Heckel et al., 2005; 
Hak et al., 2005). The CINDI campaign was performed un- 
der the auspices of the European Space Agency (ESA), of the 
International Network for Detection of Atmospheric Compo- 
sition Change (NDACC), and of the EU Framework 6 ’s Ac- 
cess to Remote Sensing of Atmospheric Constituents from 
Space - Troposat-2 (ACCENT-AT2) Network of Excellence, 
and Global Earth Observation and Monitoring of the Atmo- 
sphere (GEOMon) Integrated Project. ESA promotes ac- 
curacy of ground-based measurements that can be used for 
satellite data validation; NDACC promotes excellence in 
measurements of atmospheric composition; and GEOMon 
has been responsible for maintaining and developing net- 
works of ground-based remote sensors, in support of the 
preparation of the Global Monitoring for Environment and 
Security (GMES) Atmospheric Service. 

The CINDI campaign was organized around two main 
objectives. The first one was to determine the accuracy 
of state-of-the-art ground-based atmospheric NO 2 measure- 
ment techniques, to test these techniques under different at- 
mospheric conditions, and to improve them, in order to pro- 
vide a local reference for the validation of global satellite 
data retrievals. This responds to the need to develop and 
maintain accurate ground-based measurement capability, in- 
cluding the detection of vertical profiles of NO 2 , expressed 
e.g. in the conclusions of the ACCENT- AT2 supported work- 
shop on “Tropospheric NO 2 measured by satellites” (Piters, 
2007). The second objective was to investigate the potential 
of the different ground-based instruments for satellite data 
validation, focusing on the possibility of mapping the three- 
dimensional NO 2 distribution. The aim of CINDI is to pro- 
vide an estimate of the accuracy of retrieved tropospheric 
NO 2 columns and profiles for the participating instruments 
under different atmospheric conditions (pollution levels, me- 
teorological state) and viewing geometries and to provide a 
set of recommendations concerning instrumentation, calibra- 
tion, and measurement and retrieval strategies for the vali- 
dation of satellite-derived tropospheric NO 2 data. A better 
understanding of the variability of NO 2 could very well im- 
prove current state-of-the-art air quality models. 

These main objectives resulted in a two-phase approach. 
The first phase of CINDI, running from 8 June to 8 July, 
consisted of an extensive intercomparison campaign for 
~22 multi-axis and zenith-sky spectrometers that were lined 
up together at Cabauw. The result of this study, showing that 
the slant columns of NO 2 and O 4 measured by these instru- 
ments, when retrieved with standardized setting, are mostly 
consistent within 5 to 10%, has been published in Roscoe 
et al. (2010). In the second phase of CINDI, running from 
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9 July to 26 July, the observation modes and locations of 
the instruments were optimized for investigating the influ- 
ence of different viewing geometries and for investigating 
spatial variability on the scale of a typical satellite ground 
pixel footprint over the measurement site. 

Measuring the vertical distribution of NO 2 was of key 
importance in the CINDI campaign as its uncertainty is a 
major source of error in the retrieval of tropospheric NO 2 
columns from both ground-based and satellite based obser- 
vations. Next to the passive remote sensors like the Multi- 
AXis Differential Optical Absorption Spectrometers (MAX- 
DOAS, Honninger et al., 2004), an active remote sensing 
system consisting of a lidar instrument developed at RIVM 
(Volten et al., 2009) measured the vertical distribution of 
NO 2 in the boundary layer. Additionally four in-situ chemi- 
luminescence analysers were placed at several levels in the 
Cabauw tower, providing valuable information about the 
temporal evolution of the NO 2 vertical profile in the low- 
est 200 m, and monitoring the boundary layer development. 
Simultaneous observations at 3 m altitude provided the op- 
portunity to demonstrate that there is a bias between the in- 
situ instruments working with the widely-used molybdenum 
NO 2 to NO converter, which is also sensitive to other nitrate 
containing chemical compounds such as peroxyacetyl nitrate 
(PAN; Dunlea et al., 2007), and the systems working with the 
more recent photolytic converters which do not suffer from 
this interference effect (Fehsenfeld et al., 1990). In addition 
six experimental flights of a newly developed balloon-borne 
sonde measured the NO 2 vertical profile from the ground up 
to 5 km altitude (Sluis et al., 2010). The CINDI NO 2 pro- 
file data were intercompared to evaluate the various profiling 
techniques, pointing out different strengths and weaknesses, 
and to make suggestions for future retrieval improvement 
(see Sect. 4.4). 

Another element of CINDI was the measurement of a va- 
riety of atmospheric constituents. Besides the retrieval of 
NO 2 columns and profiles simultaneous retrievals were per- 
formed of aerosol parameters (extinction profiles, aerosol 
optical thickness), ozone (O 3 ), formaldehyde (HCHO), gly- 
oxal (CHOCHO), sulphur dioxide (SO 2 ), and water vapour 
(H 2 O). Together with the permanent measurement capability 
at the Cabauw site (see Sect. 2.2), characterizing the meteo- 
rological state of the atmosphere, the CINDI campaign thus 
provided a comprehensive observational data set of atmo- 
spheric constituents in a highly polluted region of the world 
during summertime. 

The paper is organized as follows. In Sect. 2 we describe 
the campaign design, participating instruments, and the set- 
up of the experiments. In Sect. 3 we describe the meteorolog- 
ical and air pollution conditions, the data collected and the 
measurement characteristics. Early results of the campaign 
are presented in Sect. 4, and Sect. 5 discusses the potential 
of this unique data set for future studies. The conclusions are 
given in Sect. 6 . 


2 Campaign design 

2.1 Seasonal timing 

The early summer timing of CINDI was motivated by sev- 
eral reasons. Summer is the season with mo 60, 170, and 
100 for winter, spring, and autumn, respectively), increas- 
ing the probability of cloud-free skies needed for favourable 
ground-based and satellite retrievals. During summer, con- 
vection leads to well-mixed and relatively deep boundary 
layers with NO 2 concentrations that can be well-detected by 
both MAX-DOAS and satellite instruments. A disadvantage 
of having an NO 2 measurement campaign in summer, is that 
NO x emissions are generally lower than in winter, and the av- 
erage lifetime of NO 2 is shorter. In late summer the increased 
convective activity more often leads to cumulus clouds and 
even thunderstorms, obstructing remote sensing observations 
of NO 2 . In winter, when NO 2 is mostly confined within a 
shallow boundary layer, lifting to the free troposphere is sup- 
pressed. This, in combination with the low solar elevation 
angles in the winter, would lead to more challenging circum- 
stances, in particular for satellite retrievals. 

2.2 Description of the CESAR site 

The CESAR site is located at latitude 51.971° N, longitude 
4.927° E, at about 0.7 m below sea level in an extended and 
flat polder landscape as indicated in Fig. 1 (left image). The 
Cabauw site was chosen because of its unobstructed view 
close to the horizon, its large day-to-day variability in tropo- 
spheric nitrogen dioxide and aerosols enabling the sampling 
of a wide range of pollution conditions, the absence of local 
pollution sources, the 213 m tower as depicted in Fig. 1, from 
which the planetary boundary layer can be sampled at vari- 
ous altitudes, and the excellent local support. The same site 
has been used for two previous more exploratory campaigns 
focusing on the validation of satellite observations (Brinksma 
et al., 2008; Hains et al., 2010). Cabauw is a rural site, with 
only a few pollution sources nearby. However, the wider 
vicinity is densely populated, with the city of Utrecht and a 
dense highway grid within 25 km, so that the site experiences 
recurring pollution events such as from the daily morning and 
afternoon rush hours. In addition, Cabauw is influenced by 
the transport of air pollution from emission sources further 
away, as shown in the June— July average tropospheric NO 2 
column retrieved from OMI (see Fig. 2). Northerly winds 
generally carry relatively clean air from the sea, but winds 
from any other direction are likely to result in the sampling of 
polluted air. For winds from the west to south-west, Cabauw 
is downwind of Rotterdam (40 km), Europe’s largest harbour 
and location of petrochemical plants. Inflow from the south 
to south-east carries pollution from the southern parts of the 
Netherlands, Belgium, and the industrialized German Ruhr 
area (140 to 190 km). 
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Fig. 1. Left panel: the 213 m tall meteorological tower at the CESAR site, Cabauw, The Netherlands (source: Wikipedia, http://nl. wikipedia. 
org/wiki/Meetmast_Cabauw). Right panel: aerial view of the CESAR site at Cabauw depicting the grassland polder landscape, a nearby 
village to the North, and a local secondary road to the South. Here “RSS” indicates the remote sensing site and “WP” the wind profiler site. 
Image orientation is north-south hence nearby “Zijdeweg” runs at 23° away from due north. Depicted area measures 900 x 1 100 m 2 (source: 
Google Maps, http://maps.google.nl). 


OMI tropospheric N0 2 (vl.02) 8 June - 25 July 2009 



^ __ ! 

0.00 1.50 3.00 4.50 6.00 7.50 9.00 10.5 12.0 13.5 15.0 


Fig. 2. Average tropospheric vertical NO? column during 
CINDI (between 8 June and 25 July 2009) as measured by OMI 
(DOMINO vl.02). The campaign site Cabauw is located almost 
in the middle of the Netherlands. Nearby pollution sources are the 
harbours of Rotterdam (40 km) and Antwerp (Belgium, 90 km) and 
the industrial Ruhr area in Germany (140 to 190 km). 


The three measurement locations at the CESAR site, used 
for CINDI, are shown in Fig. 1 (right image). At the remote 
sensing site (RSS) most of the MAX-DOAS and direct-sun 
instruments and the lidar systems were deployed. The nearby 
wind profiler site (WP) provided extra room for several 


MAX-DOAS instruments and additional zenith-sky viewing 
systems. The tower housed all the in-situ sensors taking in 
air at ground and elevated levels, and had MAX-DOAS and 
zenith-sky instruments and in-situ systems installed at vari- 
ous altitudes. 

2.3 Participating instruments 

The principal instruments that participated in CINDI, their 
location on site and relevant specifications are summarized 
in Tables 1 and 2. Detailed instrument descriptions are given 
in Appendix A. 

More than 20 optical spectrometers originating from 14 in- 
ternational research institutes took part in the field campaign 
and measured spectra of scattered sunlight from the ground 
for 4 to 7 weeks. From these spectra, slant columns of ni- 
trogen dioxide and several other species (see Sect. 3.3) were 
derived by means of the differential optical absorption spec- 
troscopy technique (DOAS, Platt, 1994). Slant columns of 
NO 2 and the oxygen dimer (O 2 -O 2 ) were compared in a 
formal semi-blind comparison, following NDACC criteria 
(Roscoe et ah, 2010). Most of the instruments were MAX- 
DOAS systems, capable of performing scans in elevation and 
sometimes azimuth while measuring scattered sunlight (see 
for example Honninger et ah, 2004; Wittrock et ah, 2004; 
Clemer et ah, 2010). One of the National Aeronautics and 
Space Administration (NASA) instruments was direct-sun 
viewing only, and the French National Center for Scientific 
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Fig. 3. Snapshots of the various optical instruments installed at the remote sensing site indicating the large variety of optical detection system 
designs actively involved in the campaign effort. Some images display complete stand-alone systems (e.g. WSU), others are only displaying 
the telescopes (e.g. NASA) that are connected to the spectrographs and detection systems by optical quartz fibres. 


Research (CNRS) Systeme d’Analyse par Observation Zen- 
itale (SAOZ) systems were zenith-viewing only. An illus- 
trative collection of images is shown in Fig. 3, highlighting 
the tremendous variety in the design of the optical detection 
systems. 

Most optical spectrometers deployed during CINDI con- 
sist of three basic components; the entrance optics, the spec- 
trometer and the detector. Figure 4 summarises the lay- 
out and viewing geometry of these instruments. The en- 
trance optics includes a telescope that defines the Field of 
View (FOV). It is often mounted on a motorized mechan- 
ical tracker. The telescope is connected to the spectrom- 
eter by means of quartz fibre bundles. The spectrometer 
subsequently separates the incoming light as a function of 
wavelength, and projects it on the detector. The spectrome- 
ter and detector are often housed in the same unit. Excep- 
tions to this description are the WSU, University of Toronto, 
and NIWA systems, and the six “mini-MAX-DOAS” instru- 
ments (commercially available all-in-one systems). The en- 
tire WSU system is mounted on a sun-tracker, the stationary 


NIWA system contains an optical directional scanner inter- 
nally, and the stationary system of the University of Toronto 
can switch between zenith and MAX-DOAS functionality by 
having an optical scanner fitted to its zenith entrance port. 
The telescope of the mini-MAX-DOAS systems is attached 
to the instrument and the complete instrument is mounted on 
a stepper-motor. All systems are different in terms of FOV, 
spectral range and resolution, see Table 1 . 

The RIVM nitrogen dioxide lidar with elevation scanning 
and operating from a mobile truck provided vertical pro- 
files of nitrogen dioxide at moderate resolution (Volten et al., 
2009). The Raman lidar Caeli provided vertical profiles of 
aerosol optical properties, backscatter at three wavelengths 
and extinction at two wavelengths (Apituley et al., 2009). 
Additional in-situ nitrogen dioxide instruments equipped 
with photolytic converters (Kley and McFarland, 1980) were 
operated on the ground and on the various tower platforms. 
A novel balloon-borne nitrogen dioxide sonde developed by 
KNMI (Sluis et al., 2010) was flown for the first time, provid- 
ing in-situ measurements of NO 2 at high vertical resolution. 
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Fig. 4. Left panel: MAX-DOAS instrument layout. The dashed box indicates that the three components (optics, spectrometer, detector) can 
be separate units or included in the same housing. Right panel: MAX-DOAS viewing geometry. The participating MAX-DOAS instruments 
measured at different elevation angles consecutively; tropospheric information can be retrieved using the fact that the light path through the 
stratosphere is almost constant as a function of elevation angle. The scanning is performed by rotating either a mirror or prism, or the optical 
head or telescope, or the whole instrument. 


Ozone was measured with an ozone lidar (Hoexum et al., 
2009; Apituley et al., 2010) located in Bilthoven (30 km from 
the site) and by ozone sondes launched from the CESAR site. 

TNO continuously performs aerosol observations in the 
basement of the tower where air is sampled via a common 
inlet at 60 m. A multi-angle absorption photometer (MAAP) 
was used to quantify the aerosol absorption coefficient. An 
integrating nephelometer was used for the (back-) scattering 
coefficient. To increase comparability between observations 
in (global) aerosol networks (WMO/GAW guidelines, 2003) 
prescribe that sampled aerosol is dried to relative humidities 
below 40%. Aerosol optical properties, most notably the 
scattering coefficient, strongly increase with increasing rel- 
ative humidity; thus drying frustrates comparison to aerosol 
optical properties measured at ambient conditions, e.g. re- 
motely sensed aerosol properties. To overcome this issue, 
Paul Scherrer Institute (PSI) operated a humidified neph- 
elometer that measures the aerosol (back-)scattering coeffi- 
cient at defined relative humidities between 20 and 95 %. 

The placement of the instruments over the campaign site 
was driven by various requirements. The most important 
of which was the accessibility of all desired viewing angles 
while avoiding mutual visual obstructions by the telescopes. 
Many receivers were deployed on the roof tops of four con- 
tainers or on the concrete surface next to them with most of 
the spectrometers inside. The roof-top deployment provided 
an extra 3.5 m of altitude hence overlooking most of the low 
vegetation. The second requirement was the need for tak- 
ing observations at higher altitudes which were performed 
from the tower platforms by both in-situ and remote sensing 
systems. 


In order to meet the power needs for all participating in- 
struments, an additional electrical power generator was in- 
stalled at the site. We chose to install a 100 kW diesel-driven 
generator at a fair distance from the telescopes (30 m) and 
in such an orientation that its exhaust gases were transported 
away from the site by the most frequently occurring wind 
directions (south-west to west). The nominal load on the 
generator was around 20 % of its maximum output allow- 
ing the brand-new diesel engine to run smoothly and effi- 
ciently. During the campaign no sign of contamination of 
the in-situ and remote sensing observations by the generator 
exhaust gases has been reported. 

2.4 Ancillary observations 

In support of the systems temporarily deployed at the site 
for the campaign purpose as described in Sect. 2.3, the CE- 
SAR site at Cabauw is equipped with a large suite of instru- 
ments for continuous monitoring and characterization of the 
state of the atmosphere. Standardized systems record the sur- 
face (2 m) level temperature, pressure, visibility, wind direc- 
tion and speed, and humidity. The KNMI-tower at Cabauw 
was originally built for meteorological research to establish 
relations between the state of the planetary boundary layer 
(PBL), the land surface conditions and the general weather 
situation for all seasons. It also supports aviation at Schiphol 
Airport. Hence the tower has wind, pressure and tempera- 
ture instruments installed at various levels. Nitrogen diox- 
ide and aerosols are routinely sampled close to the base of 
the tower by in-situ instruments as part of the national air 
quality network (LML) and the CESAR partnership. The 
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Table 1. Overview of CINDI participating spectrometers. “MMD” means Mini-MAX-DOAS (Appendix A). Columns denote: institute, 
field of view, spectral range and spectral resolution, coupling between telescope and spectrometer (F: multimode fibre, D: direct coupling, 
L: liquid light guide), type of detector, detector temperature and site location (RSS: remote sensing site, Twr: tower, WP: wind profiler site, 
M: mobile). Other details are given in Appendix A. 


Institute/ instrument 

FOV(°) 

range (nm) 

res (nm) 

Coupl. 

Det. 

T (°C) 

Loc. 

BIRA-IASB D. Channel 

0.8 

300-388 

0.37 

F 

CCD 

-30 

RSS 



400-700 

0.95 

F 

CCD 

-30 


BIRA-IASB MMD 

0.6 

290-133 

0.6 

D 

LinArr 

+15 

Twr 

BIRA-IASB mobile 

2.5 

270-500 

1.15 

F 

LinArr 

n/a 

M 

IUP-Bremen D. Channel 

1.2 

315-384 

0.4 

F 

CCD 

-35 

RSS 



401-573 

0.8 

F 

CCD 

-35 


CNRS SAOZ 

9.0 

270-630 

1.1 

D 

LinArr 

n/a 

WP 

CNRS mobile 

9.0 

270-630 

1.1 

D 

LinArr 

n/a 

M 

CNRS mini-SAOZ 

6.0 

270-800 

0.6 

D 

CCD 

n/a 

Twr 

GIST MMD UV 

0.5 

290-130 

0.7 

F 

LinArr 

-5 

RSS 

IUP-Heidelberg DOAS 1 

0.9 

290-790 1 

0.6 

F 

LinArr 

+10 

WP 

IUP-Heidelberg DOAS 2 

0.9 

320-165 

0.45 

F 

CCD 

-30 

WP 

INTA NEVA-II 

1.0 

387-514 

0.6 

F 

CCD 

-40 

RSS 

INTA RASAS-II 

1.0 

393-5 10 2 

0.45 

F 

CCD 

-30 

RSS 

JAMSTEC 

<1.0 

223-558 

0.7 

F 

LinArr 

+40 

RSS 

KNMI MMD UV 

0.45 

290^433 

0.6 

D 

LinArr 

n/a 

Twr 

KNMI MMD VIS 

0.4 

400-600 

0.5 

D 

LinArr 

n/a 

RSS 

Univs Leicester&Leeds 

1.0 

425-190 

0.54 

F 

CCD 

-42 

RSS 

MPIC-Mainz MMD 

1.2 

310-161 

0.6 

D 

LinArr 

-4 

Twr 

MPIC-Mainz mobile 

1.2 

320-166 

0.45-0.85 

D 

LinArr 

-4 

M 

NASA PANDORA (2, 3, 9) 

1.6 

280-530 

0.6 

F 

CCD 

+10 

RSS 

NIWA 

0.5 

389-510 

0.48 

F 

CCD 

-20 

RSS 

Univ. Toronto 

2.0 

341-545 

0.5-2. 5 

L 

CCD 

-72 

RSS 



315-380 

0.2-0. 8 





WSU MF-DOAS 

1.4 

282^198 

0.83 

D 

CCD 

-70 

RSS 


Notes: 1 IUP Heidelberg connects three MAX-DOAS spectrometers to one fibre input to cover this large wavelength range; “after 6 July: 3 25 — 44fi nm. 


wind profiling radar determines winds throughout the lower 
troposphere. A Vaisala LD-40 Ceilometer operates contin- 
uously to determine the cloud-base altitude and the PBL 
height. RIVM operated a backscatter lidar providing infor- 
mation on cloud structures and aerosol layers overhead up 
to 15 km at 5 min resolution. Finally, the total sky imager 
(TSI) operates continuously at the site generating hemispher- 
ical images of the cloud structures overhead at a 1 min inter- 
val that are analysed for total cloud cover. The CESAR site 
at Cabauw hosts a certified BSRN (Baseline Surface Radia- 
tion Network) and AERONET (AErosol RObotic NETwork) 
site that completes the suite of instruments and observa- 
tions to characterize the state of the atmosphere. The BSRN 
instruments monitor direct, diffuse, global and downward 
short-wave and long-wave radiative components (Ohmura et 
al., 1998). AERONET is a network of ground-based sun 


photometers which measure atmospheric aerosol properties 
(Holben et al., 1998). 

2.5 Set-up of the experiments 

2.5.1 Viewing in the same direction 

In the first month of the campaign, between 8 June and 
8 July, all MAX-DOAS spectrometers and the NO 2 lidar 
scanned at different elevation angles along the same west- 
north-westerly direction, around 287°. In this direction all 
instruments had an unobstructed view down to 0.5° elevation 
angle. For the instruments at the RSS this azimuth direc- 
tion was determined by a fixed point on the horizon at about 
1.2 km distance. The spread in the locations of the instru- 
ments at the RSS results in a spread in the viewing directions 
of about 1°, which is within the FOV of most of the systems. 
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Table 2. Other CINDI participating instruments, measuring NO 2 , ozone, and aerosol. Site locations are RSS: remote sensing site, Twr: tower, 
WP: wind profiler site. Other details are given in Appendix A. 


Institute 

Instrument 

Measured parameter 

Location 

RIVM 

NO2 lidar 

NO2 profile 

RSS 

RIVM 

NO x analysers with photolytic converter 

NO, NO2 concentration 

Twr (3 m/100 m), Twr (200 m) 

RIVM 

Tropospheric O3 lidar (DIAL) 

O3 profile 

Bilthoven (~30 km NE) 

RIVM 

Caeli lidar 

Profiles of aerosol, water vapour and clouds 

RSS 

RIVM 

Backscatter lidar 

Backscatter profile 

RSS 

KNMI 

NO2 sondes 

NO2 profile 

Twr 

KNMI 

O3 sondes 

O3 profile 

Twr 

IUP Bremen 

NOx analyser with photolytic converter 

NO, NO2 concentration 

Twr (3 m) 

EMPA 

NO x analyser with photolytic converter 

NO, NO2 concentration 

Twr (3 m) 

RIVM/LML 

NO x analyser with molybdenum converter 

NO, NO2 concentrations 

Twr (3 m) 

PSI 

Humidified nephelometer (WetNeph) 

aerosol (back-)scattering coefficient (RH-dependence) 

Twr 

TNO 

SMPS/APS 

particle size distribution 0.01-10 pm 

Twr 

TNO 

MAAP 

absorption coefficient 

Twr 

TNO 

Nephelometer 

aerosol (back-)scattering coefficient 

Twr 


The instruments at the tower scanned at a fixed azimuth di- 
rection around 286°, so that they viewed approximately the 
same air mass as the systems at the RSS. The MAX-DOAS 
instruments had a minimum set of prescribed elevation an- 
gles, and each vertical scan had to be obtained within 30 min. 
All instruments viewed in zenith direction in a 30-min win- 
dow around noon. Calibration of the zero degree elevation 
angle was the subject of particular attention in the early phase 
of the campaign. Roscoe et al. (2010) pointed out that this 
calibration is crucial for the interpretation of the data, espe- 
cially for the lower elevation angles, since the air mass fac- 
tors at low elevation angles change considerably with small 
changes in elevation angle. 

2.5.2 Viewing in different directions and at different 
locations 

In the second part of the campaign, 9 July to 26 July, the 
viewing directions of some of the MAX-DOAS instruments 
and the NO 2 lidar were changed in order to study (i) the con- 
sistency of air mass factor calculations and (ii) the spatial 
variability of NO 2 . The azimuth positions for this period 
are listed in Table 3. The MAX-DOAS instruments from 
BIRA-IASB and IUP Bremen and the NO 2 lidar from RIVM 
scanned several azimuthal directions quasi-simultaneously. 
In addition the three mobile systems were mounted on a car 
and driven around Cabauw on several tours. One of these 
tours, driven by the MPIC Mainz car, is shown in Fig. 5. 
This figure also indicates the large variability in tropospheric 
NO 2 columns along the route. On 15 and 16 July the NO 2 
lidar truck was placed at 2.1 km distance viewing into the di- 
rection of the tower towards the three in-situ instruments in 
order to verify the NO 2 concentration measured by the lidar. 



Fig. 5. Mobile measurements around Cabauw with the mini-MAX- 
DOAS instrument of MPIC Mainz in the morning of 12 June. 
The wind was coming from west to north-west, with a speed of 
4.5 ms -1 . The track displayed runs (anticlockwise) over major 
highways. The colours indicate measured tropospheric vertical 
NO 2 columns. 


2.5.3 Observing from different altitudes 

The zenith-viewing mini-SAOZ instrument (Sect. 2.3) was 
located at the 213 m altitude level during the entire campaign. 
When combining the measurements with the SAOZ instru- 
ment at the ground, the partial column between 0 and 213 m 
can be derived. 

During the first part of the campaign, three of the in-situ 
monitors with photolytic converters (Sect. 2.3) were placed 
close by the regular on-site in-situ monitor with molybde- 
num converter, in order to be able to investigate the biases 
between these instruments. One of the RIVM monitors with 
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Table 3. Instruments participating in the experiment with multiple azimuthal viewing directions. 


Instrument 

Viewing directions 

Period 

IUP-Bremen VIS 

17°, 62°, 152°, 197°, 242°, 287°, 332° 

12 June-21 July 

BIRA-IASB VIS 

45°, 90°, 135°, 180°, 225°, 270°, 315°, 358° 

10-22 July 

NASA P2 

66° 

14-20 July 

NASA P3 

187° 

14-20 July 

KNMI VIS 

253° 

10-31 July 

INTA NEVA 

287° 

28 June-21 July 

Univ. Toronto 

287° 

10-17 July 

JAMSTEC 

287° 

9 June-24 July 

Univs. of Leicester & Leeds 

326° 

10-22 July 

NO 2 lidar 

224°, 233°, 242°, 251°, 260°, 269°, 278°, 287°, 288°, 296°, 305°, 314° 

9-22 July 
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Fig. 6. Overview of meteorological conditions during the CINDI campaign period. From top to bottom: hours of sunshine, hours of 
precipitation, daily maximum temperature, wind direction (pointing up means winds coming from the south), and average daytime surface 
level NO 2 derived from the in-situ instruments of RIVM, IUP Bremen and EMPA. Category A days (see text) are coloured green, category B 
days are coloured yellow. 


photolytic converter was placed at 200 m during the entire 
campaign. In the second part of the campaign the other 
RIVM monitor with photolytic converter was placed at the 
100 m level. 

The three mini-MAX-DOAS instruments at the tower 
were at 20 m during the first part of the campaign, and later 
displaced at different altitude levels in order to study the con- 
sistency of the column and profile retrievals for different alti- 
tudes. A fourth mini-MAX-DOAS instrument was placed at 
20 m in the tower, viewing in the opposite direction. 


3 Campaign execution 

3.1 Meteorological conditions during CINDI 

June and July 2009 were relatively warm and sunny in the 
Netherlands. These months were 0.5 °C warmer and had 
24% more sun hours than the long-term average (1971- 
2000) at Cabauw. The number of sun hours, maximum sur- 
face temperature, precipitation and wind direction for the 
campaign period are shown in Fig. 6. More details are given 
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in Appendix B. During CINDI, there were 9 days with 10 or 
more sunny hours, where a sunny hour is defined as an hour 
with less than 50 % cloud cover and with more than 0.5 h sun- 
shine: 12, 13, 23, 24, and 30 June, and 2, 4, 16, and 25 July. 
These days were most suitable for retrievals from ground- 
based and satellite remote sensing observations over Cabauw. 
During these “category A” days (green in Fig. 6 ), the wind 
was more frequently coming from northerly directions, com- 
pared to prevailing south-westerly winds during the rest of 
the campaign. Thirteen days (yellow in Fig. 6 ) had be- 
tween 4 and 9 sunny hours. During these “category B” days 
the wind was most frequently coming from south-westerly 
directions. 

3.2 Air pollution during CINDI 

Figure 2 shows average tropospheric NO 2 columns observed 
from OMI over north-western Europe during the CINDI 
campaign (based on Boersma et al., 2007). We see that 
Cabauw is adjacent to strong sources of pollution in the 
west (Rotterdam), the south (industrial Flanders), and the 
south-east (German Ruhr Area), but that regions north of 
Cabauw are relatively clean. Surface NO 2 concentrations 
at Cabauw measured with in-situ techniques show consider- 
able day-to-day variation during CINDI, with generally low- 
est NO 2 concentrations (2 to 4ppbv) on days with precipi- 
tation, cloudiness, or prevailing northerly winds (see Figs. 6 
and 7). On “category A” days, with sunny conditions, ob- 
served surface NO 2 concentrations range from 4 to lOppbv. 
Figure 8 illustrates that the variability in NO 2 levels on 
sunny days over Cabauw is mostly driven by wind direction. 
On days with winds from southerly and westerly directions 
(12 June, 17 June, 3 July, 14-16 July), both OMI retrievals 
and CHIMERE simulations over Cabauw show tropospheric 
NO 2 columns that are approximately 2 times higher than on 
days with winds from the clean north-east (24 June, 25 June, 
30 June, 1 July). Particularly 3 and 14 July stand out as 
strongly polluted days at Cabauw, with high concentrations 
of NO 2 (more than lOppbv on 14 July when winds were 
from the polluted Ruhr area in Germany) and HCHO, and 
secondary pollutants ozone and aerosol (see Fig. 5 in Irie et 
al., 2011). Enhanced SO 2 concentrations (more than 5ppbv, 
not shown) were sampled on 5, 6 and 14 July (Irie et al., 
2011), when the air over Cabauw was predominantly from 
the nearby North Sea and northern Atlantic. Over these re- 
gions, shipping traffic is dense, suggesting that the enhanced 
SO 2 concentrations over Cabauw originated from shipping 
activity. 

Enhanced trace gas concentrations likely originating from 
biogenic emissions were also observed during CINDI. 
Episodes with enhanced HCHO concentrations were ob- 
served on several days. For example, on 25 and 26 June, 
and on 3 July 2009 HCHO concentrations inferred from 
MAX-DOAS measurements reached 4-8ppbv. On these 
days, maximum afternoon temperatures at Cabauw were 



Fig. 7. Average surface NO 2 mixing ratio during CINDI as a 
function of wind direction. Included are data between 06:00 and 
1 8:00 UTC and with wind speed exceeding 1 .5 m s _ 1 . 


between 25 and 30 °C. The strong correlation between 
HCHO concentrations and ambient temperature at Cabauw 
reported by Irie et al. (2011) suggests that the enhanced 
HCHO concentrations observed during CINDI originated 
from nearby isoprene emissions. Further discussion of the 
enhanced VOC concentrations, and their relations to sources, 
is given in Irie et al. (2011). 

3.3 Measurement characteristics 

The following parameters were inferred from the 
measurements: 

- in-situ concentrations of NO 2 at different altitudes along 
the tower; 

- slant and vertical columns of tropospheric NO 2 with the 
MAX-DOAS, direct-sun and zenith-sky instruments; 

- vertical profiles of NO 2 with a lidar and a sonde, 
and with the MAX-DOAS as inferred from the slant 
columns at different elevation angles; 

- slant columns of other species, like O 4 , HCHO, SO 2 , 
O 3 , H 2 O, and CHOCHO, with some of the MAX- 
DOAS and direct-sun instruments - for some of the 
trace gases (e.g. HCHO) also vertical columns and pro- 
file information have been derived; 

- aerosol optical depth with the MAX-DOAS instruments 
and the AERONET instrument; 

- aerosol extinction profiles from Caeli (from first princi- 
ples), and with the MAX-DOAS (as inferred from the 
O 4 slant columns); 
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Fig. 8. Tropospheric vertical NO 2 columns observed by OMI and simulated by the CHIMERE chemistry transport model on mostly sunny 
days at Cabauw during CINDI in June (left columns) and July (right columns) 2009. Grey areas in the OMI (DOMINO vl.02) satellite 
observations generally indicate cloud radiance fractions >0.5. The white trajectory shows the path that an air parcel has travelled from 
03:00 UTC until arriving at Cabauw at 12:00 UTC (source: http://www.nilu.no/trajectories/modeldata/2009/cabauw). The OMI observations 
are performed between 1 1:00 and 14:00 UTC, CHIMERE simulations are shown for 1 1:30 UTC. 


- aerosol extinction at the surface with in-situ (hu- 
midified) nephelometers and absorption photometers 
(MAAP) 

- ozone profiles with sondes and an ozone lidar. 

Figure 9 gives an overview of when the measurements were 
performed, detailed per instrument type. The NO 2 in-situ 
instruments from IUP Bremen, EMPA and RIVM measured 
almost continuously from 12 June to 10 July, with one in- 
strument at 200 m altitude and the rest at 3 m altitude. After 
3 July the EMPA instrument left, and from 10 July onwards 
one of the in-situ instruments was placed at the 100 m level, 
so that the concentrations at three different levels were mea- 
sured simultaneously (see Fig. 10). The in-situ instruments 


were the only instruments measuring during night and day. 
The majority of the MAX-DOAS instruments took measure- 
ments between 13 June and 21 July. Slant column densities 
were derived for the whole period with an estimated preci- 
sion of 5-10 % for most trace gases (see Roscoe et al., 2010; 
Pinardi et al., 2012). 

From the slant columns the tropospheric column amounts 
and some information on vertical profiles can be derived 
(Honninger et al., 2004; Wittrock et al., 2004). Different 
retrieval methods have been proposed and tested to convert 
slant column densities to vertical column densities and on- 
ward to vertical profiles. These methods are subject to in- 
tercomparison as briefly described in Sects. 4.3 and 4.4 and 
in more detail in Wittrock et al. (2012). All such methods 
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Fig. 9. Number of instruments which have collected data on each day during the CINDI campaign and for each of the instrument types. The 
day numbers in June and July are at the top and bottom of the figure. The yellow and green days are days of category A and B, respectively 
(see Sect. 3.1). 



Fig. 10. In-situ measurements of NCD taken at different altitude 
levels in the Cabauw tower. All measurements are performed with 
NO x analysers with photolytic converters. The upper graph shows 
the absolute values of the measurements at each level, the lower 
graph the differences. 


are sensitive to meteorological conditions, in particular the 
presence of high aerosol or cloud coverage, and therefore the 
retrieval study focused on selected sunny days. 

Zenith only measurements are available between 9 June 
and 10 July. The University of Toronto DOAS instrument 
pointed towards zenith during the first part of the campaign. 
The CNRS SAOZ instruments were always pointing towards 
zenith. Combining the surface-based SAOZ instrument with 
the mini-SAOZ on top of the tower enabled to derive partial 
columns for the first 213 m. 

Mobile measurements were performed by two MAX- 
DOAS instruments and one SAOZ instrument on 15 days 
between 12 June and 16 July, indicated in Fig. 9, on differ- 
ent tracks mostly within 30 km around Cabauw (see Fig. 5). 
Mobile measurements were also performed on longer tracks 
between Cabauw and Paris (500 km) on 6 and 10 July, be- 
tween Cabauw and Mainz (400 km) on 13 June, and between 
Brussels and Cabauw (150 km) on 13 July. The NO 2 lidar 
was operated on 28 days, indicated in Fig. 9. It measured 
NO 2 profiles up to 2.5 km (sometimes up to 4 km) the max- 
imum altitude being determined from the loss of return sig- 
nal at higher altitudes, depending on atmospheric conditions 
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like the presence of aerosols, water vapour, or clouds. The 
measurements were performed all day. Only the presence of 
precipitation or dense fog prohibited measurements. A verti- 
cal NO 2 profile is obtained by scanning a range of elevation 
angles, each elevation angle being sensitive to a different al- 
titude. The lidar measured NO 2 concentrations over a ver- 
tical range that depended on altitude, from 12 m range for a 
measurement at 12 m altitude to 1800 m for a measurement 
at 2200 m altitude. These concentrations are very accurate 
when compared to the in-situ monitors at 0 , 100 and 200 m 
altitude (Berkhout et al., 2012). 

An experimental NO 2 sonde was launched on 6 days, in- 
dicated in Fig. 9, at 10:30 UTC (Sluis et al., 2010). These 
sondes delivered high-resolution in-situ NO 2 profiles up to 
5 km. On the same 6 days, at 13:30 UTC, ozone sondes were 
launched. Both ozone and NO 2 sondes had a radiosonde 
attached to them, recording temperature, pressure, and rel- 
ative humidity. In Bilthoven, 25 km north-east of Cabauw, 
an ozone lidar measured tropospheric ozone profiles on the 
same 6 days, every 20 to 30 min. 

Aerosol properties were measured by a suite of instru- 
ments. Measurements of aerosol optical depth (AERONET, 
only during daylight hours) and scattering and absorption co- 
efficients were performed by in-situ systems during the entire 
duration of the campaign. The in-situ particle number size 
distribution for the smallest particles (0.01 to 0.5 pm), the 
SMPS, was not operational between 14 June and 5 July and 
after 19 July. 


4 Early results 

The extensive data set obtained during CINDI enables us to 
perform in-depth intercomparisons to assess the accuracy of 
the NO 2 measuring instruments. These studies will point to 
the possible origins of differences in the derived geophysi- 
cal parameters, and whether they are related to instrumental, 
algorithm, or interpretation differences. 

4.1 Accuracy of slant columns of NO 2 , O 4 , and HCHO 

The first two weeks of the campaign, a semi-blind intercom- 
parison was performed involving the 22 MAX-DOAS and 
zenith-sky instruments (Roscoe et al., 2010). These instru- 
ments pointed in the same direction and scanned almost the 
same air mass every 20-30 min. The wavelength range and 
other settings describing the DOAS analysis to obtain slant 
columns of NO 2 and O 4 were prescribed to minimise differ- 
ences caused by these algorithms. Although no systematic 
investigation was performed to evaluate the impact of not 
following these recommendations, it is known from past ex- 
perience (see e.g. Vandaele et al., 2005; Roscoe et al., 1999) 
that the use of different NO 2 cross-sections can be a signif- 
icant source of divergence. E.g. retrievals performed with- 
out consideration of the NO 2 temperature dependence may 
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Fig. 11. Straight-line slopes and their standard errors and inter- 
cepts of NO 2 slant columns for each instrument participating in the 
semi-blind intercomparison at visible wavelengths against those of 
a reference data set. Colours refer to elevation angles shown in 
the top right legend. Note that the MPIC Mainz instrument used a 
non-standard wavelength range for spectral analysis because of the 
limited range of the instrument (Fig. 6 from Roscoe et al., 2010). 


introduce a bias as high as 20 % on the slant columns. Com- 
paratively the sensitivity to the wavelength range used for 
NO 2 fitting is smaller, although bias can still be produced due 
to possible enhanced interference with water vapour absorp- 
tion, the Ring effect or the wavelength dependency of the air 
mass factor. The 30-min averages of the slant column densi- 
ties of NO 2 and O 4 matched within 5 % for most of the par- 
ticipating instruments and all elevation angles (see Fig. 11). 
Restricting the comparison to zenith-sky data, it was con- 
cluded that all instruments meet the criteria for endorse- 
ment by NDACC (http://www.ndsc.ncep.noaa.gov/organize/ 
protocols). From the study, it became clear that a proper cal- 
ibration of the elevation angles is crucial for the interpreta- 
tion of the MAX-DOAS measurements. Also, the temporal 
variability of NO 2 was quite large over typical measurement 
sampling time scales (6 to 30 min, see Table Al). There- 
fore a significant part of the scatter in the comparisons was 
caused by time differences in the measurements. A better 
synchronisation of the measurements is necessary for future 
intercomparison exercises. 

Slant columns of HCHO have been compared for a number 
of MAX-DOAS instruments. In order to reduce the retrieval 
impact, common retrieval settings (i.e. wavelength range, 


www.atnios-meas-tech.net/5/457/2012/ 


Atmos. Meas. Tech., 5, 457-485, 2012 


470 


A. J. M. Piters et al.: The CINDI campaign: design, execution and early results 


HCHO DSCD at 4° elevation on 2 July 2009 


80 ■ 


cn—70 
E 


60 


50 



• BIRA 





• IUP Bremen 





• IUP Heidelberg 





• JAMSTEC 










• WSU 



• • 


TORONTO 



• 


MPI Mainz 



• • 


• NIWA 


• 

a J 


40 ■ 


30 


• i 
i • 


1 ! * : ! M : ! J I i 


• f • i * 

: s : : I 


: • 


10 


15 


20 



Fig. 12. Formaldehyde (HCHO) differential slant column density intercomparison for the data obtained on 2 July 2009 by all systems listed 
in the legend. The left panel shows the evolution over the day of the 30-min averaged slant columns retrieved from the 4° elevation angle 
data, while the right panel shows their scatter plot with respect to the Bremen instrument (Pinardi et al., 2012). 
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Fig. 13. Top panel: comparison of four in-situ atmospheric NO 2 sampling instruments. The green line represents the instrument with a 
conventional hot-plate molybdenum converter, whereas the black, red and blue lines represent the instruments with a photolytic converter, 
from RIVM, Empa, and IUP Bremen, respectively. The middle and lower panels show the absolute and relative difference (enhancement) 
between the mixing ratio derived from the instrument with molybdenum converter and the average mixing ratio derived from the three other 
instruments. 


cross-sections and DOAS settings) were used by the different 
groups and 30-min averages were taken to reduce differences 
in temporal sampling. An example of the comparison for 
2 July 2009 and one elevation angle, including 8 different in- 
struments, is presented in Fig. 12. Over the whole campaign, 
the HCHO slant columns retrieved by the different instru- 
ments from the off-axis elevations generally agreed within 


15 % (Pinardi et al., 2012). This is a very encouraging result 
considering that this represents a first attempt to intercom- 
pare HCHO from MAX-DOAS instruments. 
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Fig. 14. Tropospheric NCb columns derived from direct-sun (red) and MAX-DOAS (blue) measurements from the BIRA-IASB and 
WSU MAX-DOAS instruments, and from zenith-sky measurements (black) from the CNRS/SAOZ instrument, on 30 June and 4 July 2009. 
The measurements from BIRA-IASB and WSU were averaged for data points within 10 min of each other. 


4.2 Comparing NO 2 concentration from in-situ 
monitors and lidar 

Five chemiluminescence in-situ monitors, one with a molyb- 
denum converter from RIVM/LML, and four with a pho- 
tolytic converter, from EMPA, RIVM (2 instruments), and 
IUP Bremen participated in the campaign. Figure 13 shows 
the measurements of four of these instruments over five days 
during the campaign. The instrument with molybdenum con- 
verter clearly has a positive bias with respect to the instru- 
ments with photolytic converters. This is most probably 
caused by the well-known interference with other species 
(e.g. with PAN and FINO3, see Winer et al., 1974). The bias, 
which is up to 5 ppbv, occurs both during day and night, but 
in relative terms is largest during the day, when solar UV 
induced atmospheric chemistry generates most of the inter- 
fering species. 

On 15 and 16 July the RIVM NO2 lidar was placed at 
a distance of 2. 1 km from the tower and pointed at the in- 
situ instruments placed at three different altitudes. The hor- 
izontal distance of the NO2 lidar to the tower was chosen 
such that the air mass measured by the lidar contained the air 
mass sampled by the in-situ monitors. Berkhout et al. (2012) 
show that the lidar measurements follow the in-situ measure- 
ments within 5% and capture variations down to a 10 min 
time scale. 

4.3 Comparing tropospheric NO 2 columns 

A first comparison was performed between the tropospheric 
vertical columns from MAX-DOAS, zenith-sky and direct- 
sun measurements (Spinei et al., 2012), see for example 
Fig. 14. The slant columns were converted to vertical 
columns, using air mass factors calculated with radiative 
transfer models (RTMs). The influence from stratospheric 
NO2 is negligible for MAX-DOAS measurements, but it has 
to be taken into account when tropospheric NO2 columns are 


derived from zenith-sky or direct-sun total column observa- 
tions. In this case, the stratospheric NO2 column is deter- 
mined from twilight observations according to methods in 
use within NDACC (e.g. Goutail et al., 1994; Hendrick et 
al., 2004; Chen et al., 2009). A linear regression between 
the BIRA-IASB and WSU data sets was performed showing 
agreement within 5 % in all cases (direct-sun: slope 0.96, 
intercept —2.36 x 10 14 molec cm -2 , R 2 0.98; MAX-DOAS: 
slope 0.99, intercept —8.9 x 10 13 molec cm -2 , R 2 0.92). 
Therefore, the measurements from BIRA-IASB and WSU 
were averaged for data points within 10 min of each other. 
SAOZ tropospheric columns were calculated using air 
mass factors based on MAX-DOAS derived profiles. The 
tropospheric columns of the three types of measurements 
agree on average within 30%. Vlemmix et al. (2011) com- 
pared partial tropospheric vertical NO2 columns from the 
MAX-DOAS with those derived from the NO2 lidar (up 
to the maximum altitude of the lidar measurement). They 
found a correlation of 0.78 and an average difference of 
0.1 x 10 15 molec cm -2 . 

4.4 Comparing vertical profiles of NO 2 

Several different methods have been developed to retrieve 
vertical profiles from MAX-DOAS measurements (see, e.g. 
Wittrock et al., 2004; Heckel et al., 2005; FrieB et al., 2006; 
Irie et al., 2008, 2011; Clemer et al., 2010; Li et al., 2010; 
Wagner et al., 2011; Vlemmix et ah, 2011). The strength and 
weaknesses of these methods are still not well understood 
and one of the objectives of CINDI was to further develop 
and improve ground-based profiling with MAX-DOAS. The 
number of degrees of freedom for MAX-DOAS profile re- 
trievals is between 2 and 5, depending on clouds and aerosols 
and on the profile shape itself (see, e.g. Clemer et ah, 2010; 
Vlemmix et ah, 2011). The different retrieval methods were 
applied to well-prepared data sets of NO2 slant columns and 
compared to each other (Wittrock et ah, 2012). The first data 
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Fig. 15. N02-profiles for 25 June 2009, 10:50UTC for four dif- 
ferent MAX-DO AS retrieval algorithms. In green are the profiles 
derived with the optimal estimation method, in red are the profiles 
derived with a simple parameterization method. In-situ data at 3 and 
200 m levels are indicated in blue. 


set consisted of modelled NO 2 slant columns for a total of 
16 different scenarios (eight NO 2 profiles and two aerosol 
scenarios). Differences in results for this comparison can 
only be due to the retrieval techniques. The second com- 
parison is on the retrieval of real data using fixed settings for 
all optimal-estimation-like retrieval techniques. For this ex- 
ercise only data from the Bremen instrument have been anal- 
ysed to avoid differences in NO 2 results due to e.g. different 
timing of the instruments. For the comparison with comple- 
mentary measurements from in-situ instruments, NO 2 -lidar 
and NO 2 sondes, “best settings” for each MAX-DOAS were 
chosen by the different groups by using their own slant col- 
umn results. Flere, “best settings” means that all groups us- 
ing optimal estimation (Rodgers, 2000) to retrieve the pro- 
file information usually apply different a priori information 
and also different numerical parameters to avoid physically 
meaningless results. 

The model study showed that for the most important pa- 
rameters - the tropospheric vertical column and the mixing 
ratio close to the surface - the different MAX-DOAS meth- 
ods (including the non-optimal-estimation methods) are able 
to reproduce the true numbers within 20 % even for higher 
aerosol (AOD = 0.54) and more difficult scenarios like up- 
lifted NO 2 layers. As expected, the retrieval in the visi- 
ble (around 450 nm) is more stable and results are closer 
to reality than in the UV (around 350 nm). When compar- 
ing the different retrieval algorithms using measured NO 2 
slant columns obtained by the Bremen instrument the agree- 
ment is very good (see e.g. Fig. 15). For surface concen- 
trations and tropospheric columns the correlation is usually 
better than 0.9, irrespective of atmospheric conditions, with 
slopes within 15 % of unity. Comparison to complementary 
data from the NCh-lidar and in-situ instruments show rea- 
sonable agreement for surface layer mixing ratios on sunny 
days (correlation R 2 > 0.7) with slopes within 25 % of unity. 


Fig. 16. NO 2 surface concentrations for in-situ (IUP Bremen, red), 
lidar (RIVM, blue), and MAX-DOAS (IUP Bremen, green). 


For atmospheric conditions with a very low boundary layer 
height of less than 300 m the MAX-DOAS profile retrievals 
significantly underestimate the “real” NO 2 concentration by 
up to 50 % as shown in Fig. 16 for 23 and 24 June. 

4.5 Comparing vertical profiles and optical depth of 
aerosol 

Aerosol extinction vertical profiles were retrieved from 
MAX-DOAS measurements by five groups (BIRA-IASB, 
IUP Heidelberg, JAMSTEC, KNMI and MPIC Mainz). 
While some of the retrieval algorithms retrieve full vertical 
profiles using optimal estimation, others use more simple 
least squares algorithms to retrieve one or two parameters 
(e.g. aerosol optical thickness and layer height). Despite the 
large differences in the retrieval approaches, a comparison 
with ceilometer measurements demonstrates that the general 
structure of the boundary layer from MAX-DOAS observa- 
tions of O 4 compares favourably (see Fig. 17). Reasonable 
agreement is achieved between the resulting aerosol optical 
thickness from MAX-DOAS and AERONET measurements, 
with regression coefficients better than 0.8 and slopes rang- 
ing between 0.7 and 1.2. MAX-DOAS retrievals are less 
sensitive to aerosol layers at higher altitudes, as was seen by 
a comparison with aerosol extinction profiles from Caeli. A 
publication of the results from the aerosol intercomparison is 
in preparation (FrieB et ah, 2012). 

4.6 Comparing surface aerosol extinction 

Zieger et al. (2011) compared the lowest levels of the sur- 
face aerosol extinction profiles retrieved by four MAX- 
DOAS instruments, using different algorithms, with in-situ 
surface values of the aerosol extinction. These in-situ val- 
ues were determined from the combination of measurements 
of the scattering coefficients at different relative humidi- 
ties and the dry absorption coefficient. Zieger et al. (2011) 
found good correlation between the MAX-DOAS instal- 
ments and the Caeli lidar with in-situ extinction values, but 
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Fig. 17. Aerosol profile intercomparison for 3 July (FrieB et 
al., 2012). The top plot shows the backscatter profiles from the 
Ceilometer at Cabauw. The second plot shows the same data with 
MAX-DOAS averaging kernels applied. The plots below that show 
the aerosol extinction profiles of the MAX-DOAS instruments from 
BIRA-IASB, IUP Heidelberg, JAMSTEC and MPIC Mainz (top- 
down), respectively. 


the MAX-DOAS data are generally significantly larger than 
the corresponding in-situ values by a factor ~1.5 to 3.4. The 
Caeli lidar values (extrapolated to the ground using the mea- 
sured backscatter signal) were a factor of 1.12 to 1.76 larger 
than the in-situ values. Several hypotheses were brought 
forward, for example the limited vertical resolution of the 
MAX-DOAS retrieval (especially during lofted layers) and 
possible particle losses in the in-situ inlet system. Zieger et 
al. (2011) found that the agreement is better for low AOD and 
low PBL cases (see Fig. 18). The retrieval for some of the in- 
struments also improved when ambient in-situ measurements 
of the asymmetry factor and the single scattering albedo were 
used as input. 


5 Outlook and potential of the CINDI data set 

Section 4 highlights the most important results from the 
CINDI campaign until now. More studies with the compre- 
hensive CINDI data set are anticipated. A main objective of 
the CINDI campaign is to determine the potential of the ob- 
servations by the atmospheric NO 2 measuring ground-based 
instruments for validation of tropospheric NO 2 columns de- 
rived from satellite observations by, e.g. OMI aboard EOS- 
Aura, SCIAMACHY aboard Envisat, and GOME-2 aboard 
METOP. Here the most important aspect is the accuracy of 
the determined tropospheric columns and vertical profiles, 
which is obtained by quantifying the performance of the 
measurements and retrieval methods under different atmo- 
spheric conditions. The CINDI data set provides an excel- 
lent opportunity to study this. The experiment with differ- 
ent viewing directions, as described in Sect. 2.5.2, will be 
used to study the overall consistency of air mass factor calcu- 
lations, especially related to viewing geometry, aerosol and 
cloud properties and horizontal and vertical distributions of 
trace gases. By combining measurements taken under differ- 
ent viewing directions, elevation angles, and different rela- 
tive azimuth angles, inconsistencies in the assumptions on 
aerosol and boundary layer height can be further investi- 
gated. Another valuable source of information on the consis- 
tency of air mass factors can be derived from the comparison 
between the direct-sun and the scattered light measurements. 

An important aspect to be taken into consideration when 
comparing ground-based and satellite -borne remote-sensing 
observations is their difference in sensitivity to trace gas 
abundances as a function of the altitude at which it resides. 
This can be studied by theoretical simulations using a ra- 
diative transfer model, and the CINDI data can be used to 
verify these simulations by combining the MAX-DOAS and 
zenith-sky measurements at different viewing geometries and 
under different atmospheric conditions with the in-situ mea- 
surements at different altitudes. 

Another important aspect is the difference in sampling of 
the air mass in the respective field of view. In the horizontal 
spatial domain the satellite observes ground-pixel footprints 
at the 30 to 100 km scale, while the ground-based MAX- 
DOAS horizontal field of view is at most a few km in the 
viewing direction and a few 100 m in the across viewing di- 
rection. The horizontal variability in NO 2 along a satellite 
ground pixel can be quite large, as can be seen from the mo- 
bile measurement track in Fig. 5, which spans an area of ap- 
proximately 30 x 40 km 2 and shows a range in NO 2 values 
of a factor of 6; this variability is averaged over the satellite 
ground pixel footprint. 

The mobile measurements will be combined with the mea- 
surements from Cabauw in several viewing directions and 
with air-quality models to map typical horizontal distribu- 
tions of NO 2 and quantify the effect of such distributions in 
the interpretation of satellite measurements. 
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Fig. 18. Ambient aerosol extinction coefficient retrieved from MAX-DOAS observations versus in-situ measurements brought to ambient 
conditions. The colour denotes the AOD measured by the Cimel sun photometer (AOD interpolated in accordance with the appropriate 
wavelength; grey points are times with no sun photometer measurements). The solid black line represents a bivariate linear regression 
including weights (with calculated uncertainty of slope and intercept). The dashed line is the l:l-line (Fig. 8 from Zieger et al., 2011). 


During thirty-one days the CINDI campaign experienced 
sunny conditions during one or more daily satellite over- 
passes, as detailed in Appendix C. These satellite overpasses 
can be used to study the consistency of the satellite retrievals 
and the assumptions about aerosol and NO 2 profiles with the 
results derived from CINDI data. 

This campaign should result in recommendations for op- 
timal instrumental and retrieval settings, optimal observa- 
tion strategies, and optimal intercomparison methods for the 
ground-based systems involved. 

6 Conclusions 

The CINDI campaign has been very successful in achieving 
its observational and scientific objectives. A large data set 
of continuous ground-based in-situ and remote sensing ob- 
servations of nitrogen dioxide, aerosols and other air pollu- 
tion constituents has been collected under various meteoro- 
logical conditions and under various air pollution loadings. 
The campaign included many cloud-free periods: nine days 
were labelled as “category A” days with at least 10 sunny 
hours, and 13 days as “category B” days with at least 4 sunny 


hours. The day-to-day variability of NO 2 for sunny days 
was mostly driven by wind direction, with cleaner air com- 
ing from northerly directions. The CINDI campaign expe- 
rienced periods of enhanced SO 2 with winds coming from 
the North Sea, possibly originating from ship emissions, and 
periods of enhanced HCHO with warm winds coming over 
land, suggesting an origin of biogenic emissions (Irie et al., 
2011). The skies above Cabauw were observed in multi- 
ple azimuthal and elevation directions by the MAX-DOAS 
and direct-sun instruments from morning to evening twilight. 
First detailed comparisons performed with the CINDI data 
show that 

- MAX-DOAS slant column measurements of NO 2 and 
O 4 agree within 5-10 % (Roscoe et al., 2010); 

- MAX-DOAS slant column measurements for HCHO 
agree within ~15 % (Pinardi et al., 2012); 

- the bias that in-situ NO 2 measurements derived with a 
molybdenum converter suffer from interfering species 
can be as large as 5 ppbv during day time (Sect. 4.2); 

- the NO 2 lidar measurements agree with in-situ measure- 
ments within 5 % (Berkhout et al., 2012); 
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- tropospheric columns and surface values of NO 2 from 
MAX-DOAS NO 2 profile retrievals agree within 15% 
between different instruments, and within 25 % with the 
NO 2 lidar data and in-situ NO 2 data obtained at differ- 
ent altitudes, except for the situation of very shallow 
boundary layers, where the MAX-DOAS algorithms 
underestimate the NO 2 values by up to 50% (Wittrock 
et ah, 2012); 

- MAX-DOAS aerosol optical depth is in reasonable 
agreement with the AERONET measurements, with re- 
gression coefficients better than 0.8 and slopes ranging 
between 0.7 and 1.2 (FrieB et ah, 2012, ; Sect. 4.5); 

- MAX-DOAS aerosol extinction values are generally 
well correlated with, but a factor ~ 1.5 to 3.4 larger than, 
the corresponding in-situ values at the surface (Zieger et 
ah, 2011). 

Other studies that are currently being performed include 
the comparisons of tropospheric NO 2 columns, aerosol pro- 
files, ozone profiles, and HCHO columns, the spatial vari- 
ability of NO 2 , and the application to satellite data validation. 

It is expected that the studies performed during the CINDI 
campaign will result in increased knowledge about the per- 
formance of ground-based remote sensing instruments re- 
garding the accuracy with which NO 2 and aerosol infor- 
mation in terms of vertical profiles and tropospheric/total 
columns can be derived. The intercomparisons will re- 
sult in harmonization of retrieval settings and observation 
methods, and in recommendations for building the networks 
of ground-based systems urgently needed for satellite data 
validation. 

Appendix A 

Descriptions of participating instruments 

This Appendix gives detailed descriptions of most of the par- 
ticipating instruments. Table 1 in Sect. 2.3 lists some more 
details of the spectrometers, in particular on field of view, 
wavelength range and resolution. Table Al lists the elevation 
angles, and integration times used during CINDI for each of 
the systems. Table 2 in Sect. 2.3 lists the location and mea- 
surement parameters of the other participating instruments. 

Al MAX-DOAS instruments 

BIRA-IASB MAX-DOAS 

The BIRA-IASB MAX-DOAS instrument is a dual- 
channel system, with a UV and a VIS spectrome- 
ter, an ORIEL MS260i with a grating consisting of 
1200 grooves mm -1 and an ORIEL MS 127 with a grating 
of 600 grooves mm -1 , respectively. For both the UV and 
the VIS spectrometer, the instrumental response function is 


close to a Gaussian. The output of the spectrometers is con- 
nected to two low-noise thermo-electrically cooled CCD de- 
tector systems with 2048 x 512 and 1340 x 100 pixels for 
the UV and VIS, respectively. The optical head, mounted on 
a sun-tracker, can collect direct-sun and scattered light at var- 
ious elevation (0 to 90°) and azimuth angles (0 to 360°). The 
optical head also includes a 6-position filter wheel equipped 
with transmission diffuser plates and neutral density filters. 
A full description of the instrument can be found in the pub- 
lication by Clemer et al. (2010). 

BIRA-IASB mobile DOAS 

The BIRA-IASB mobile-DOAS instrument (Merlaud et al., 
2012) is based on two similar compact Avantes spectrometers 
installed on a car. The entry slit is 50 pm, the focal length 
75 mm and the grating is a 6001mm -1 , blazed at 300 nm. 
The CCD detector is a Sony 2048 linear array with a Deep- 
UV coating for signal enhancement below 350 nm, 40 pho- 
tons per count sensitivity at 600 nm and a signal-to-noise ra- 
tio of 200. An optical head, mounted on the car window, 
holds the two telescopes with fused silica collimating lenses. 
One telescope points zenith while the other is directed 30° 
above the horizon. Two 400 pm chrome plated brass op- 
tical fibres connect the telescopes to the spectrometer. A 
GPS antenna is used for geolocation referencing of the mea- 
surements, the whole set-up is powered by the car battery 
12 V through an inverter. While measuring, the instrument 
is recording spectra continuously and simultaneously from 
the two directions. The instrument was operated on the local 
roads around Cabauw and on the highways between Utrecht 
and Rotterdam. 

BIRA-IASB, KNMI, MPIC Mainz and GIST mini 
MAX-DOAS instruments 

Six mini-MAX-DOAS instruments were present during the 
campaign. The mini-MAX-DOAS is a commercial system 
from Hoffmann Messtechnik GmbH, consisting of an Ocean 
Optics spectrograph, diffraction grating, and a CCD detec- 
tor. The whole instrument is contained in a compact metal 
box (13 x 19 x 14 cm 3 ) mounted on a stepper motor, allow- 
ing the change in viewing elevation angle. It has an active 
Peltier cooler. The six instruments differ in their spectro- 
graphs, so that the wavelength ranges, resolution and field of 
view are slightly different, see Table 1 in Sect. 2.3. 

INTA MAX-DOAS NEVA-II 

NEVA-II is a robust house-made instrument developed for 
continuous operation in continental Antarctica. It is based 
on a HAMAMATSU S703 1-1008 (256 rows x 1024 pixels) 
CCD sensor. The read-out electronics is designed and de- 
veloped at INTA, achieving a typical CCD node sensitivity 
of 2.2 pV/e - and typical read-out noise of 8e - operating at 
—40 °C. The spectrograph is a TRIAX 180 (Czerny-Turner) 
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Table Al. Elevation angles, integration times and full scan durations of the MAX-DOAS systems during CINDI. 


Institute/instrument 

Elevation angles 

IT 1 

full scan 

BIRA-IASB D. Channel 

1°, 2°, 3°, 4°, 5°, 8°, 10°, 15°, 30°, 90° 

0.1 to 60s 

1 5 min 

BIRA-IASB MMD 

-0.5°, 0°, 2°, 4°, 6°, 8°, 10°, 15°, 30°, 85° 

1 min 

15 min 

BIRA-IASB mobile 

30°, 90° 

1 min 

2 min 

IUP-Bremen D. Channel 

20 angels 

40 s 

1 5 min 

GIST MMD UV 

0°, 1°, 2°, 3°, 4°, 5°, 8°, 10°, 15°, 30°, 90° 

2 min 

20 min 

IUP-Heidelberg DOAS 1 

2°, 4°, 8°, 15°, 30°, 90° 

1 min 

6 min 

IUP-Heidelberg DOAS 2 

2°, 4°, 8°, 15°, 30°, 90° 

1 min 

6 min 

INTA NEVA-II 

5°, 9°, 16°, 31°, 90° 

0.1 s 

20 min 

INTA RASAS-II 

3°, 4°, 5°, 8°, 10°, 15°, 30°, 90° 

10s 

15 min 

JAMSTEC 

2° 2 , 4°, 8°, 15°, 30°, 90° 

40 ms /5 min 

30 min 

KNMI MMD UV 

0°, 2°, 4°, 8°, 15°, 30°, 85° 

2 min 

1 5 min 

KNMI MMD VIS 

0°, 2°, 4°, 8°, 15°, 30°, 90° 

1 min 

15 min 

Univs Leicester&Leeds 

1°, 2°, 3°, 4°, 5°, 6°, 8°, 10°, 15°, 30°, 90° 

1 to 2 min 3 

20 min 

MPIC-Mainz MMD 

-5°, -0.5°, 0°, 2°, 4°, 6°, 8°, 10°, 15°, 30°, 85° 

1 min 

1 5 min 

MPIC-Mainz mobile 

22°, 45°, 90° 

1 min 

3 min 

NASA PANDORA (2,3) 

0°, 15°, 30°, 45°, 60°, 75°, 80°, 82°, 85°, 86°, 87°, 88°, 89° 

300 ms/20 s 

5 min 

NIWA 

1°, 2°, 4°, 8°, 10°, 15°, 30°, 90° 

0.016 to 20 s/120 s 

17 min 

Univ. Toronto 

2°, 4°, 8°, 15°, 30°, 90° 

15 ms to 30 s/5 min 

30 min 

WSU MF-DOAS 

3°, 4°, 5°, 8°, 10°, 15°, 30°, 90° 

0.63 s, 1.2 s 4 

10 min 


Notes: 1 The MAX-DOAS systems add up a number of spectra at smaller integration times to avoid saturation. The integration time reported here is the total integration time for one 
elevation angle. Where two times (or time ranges) are given, the smaller one is the typical integration time (or time range) of a single spectrum; 2 3° between 8 and 21 June; -'2 min 
for each angle submitted to the comparison exercise, 1 min for all other angles; 4 single scan duration for vis and UV respectively. 


holding a 1200 grooves mm -1 holographic grating. The sys- 
tem is located in a waterproof housing stabilized at 23 °C. 
Dry Nitrogen is supplied to keep the humidity inside below 
5 %. Depolarized light reaches the spectrograph through a 
10 m fused silica fibre bundle. The optical telescope is de- 
signed to stand 1° field of view. Tilt telescope movement 
is controlled by thermally-controlled stepping motors. Mea- 
surements were performed in a continuous mode whenever 
SZA < 98°. NEVA-II was operating in zenith mode around 
twilight (SZA > 80° at dawn and SZA> 70° at dusk) and in 
off-axis mode the rest of the day. 

INTA MAX-DOAS RASAS-II 

RASAS-II is based on an Andor IDUS UV CCD cam- 
era (model DU420A-BU) and a Shamrock SR163i 
(Czerny-Turner) spectrograph holding a grating of 
600 grooves mm~ 1 . The CCD was operating at — 40 °C 
to minimize dark current of the detector. The housing 
and optics are similar to those of the previously described 
NEVA-II. The fibre telescope pan and tilt positioner is a 
Direct Perception PTV-D47. During the first period of 
CINDI the instrument operated in the visible range for NCD 
measurements. From 6 July onwards the spectral range was 
shifted to the UV for HCHO and BrO measurements. 


IUP-Bremen MAX-DOAS 

The instrument from IUP Bremen (Wittrock et al., 2004) 
consists of two grating spectrometers - one for the UV and 
one for the VIS - and a separate telescope unit connected 
to the two channels via a Y-shape quartz fibre bundle. The 
spectrometers are temperature stabilized to avoid wavelength 
drifts. The 20 m quartz fibre bundle with 76 cylindrical sin- 
gle fibres, each of which has a diameter of 150 pm, efficiently 
depolarizes the incoming light and also provides flexibility 
for the instrument set-up. The telescope housing is based on 
a commercial watertight box, where optical components like 
a lens restricting the FOV to 1.2°, a HgCd lamp for spectral 
calibrations and a commercial video camera are fitted in. The 
box is mounted on a pan-tilt-head which is able to turn in al- 
most any direction with high speed. In contrast to most of 
the other instruments, scans in different azimuth angles were 
carried out during the campaign. The power consumption 
for the whole set-up is about 1000 VA peak, space required 
inside 1 m 2 . 

JAMSTEC MAX-DOAS 

The JAMSTEC MAX-DOAS instrument consists of a tele- 
scope unit placed outdoors and a spectrometer unit in- 
doors. These were connected by a multimode step-index 
fibre (NA = 0.22, pure silica). A miniaturized UV/visible 
spectrometer (Ocean Optics USB4000) equipped with a 


Atmos. Meas. Tech., 5, 457-485, 2012 


www.atmos-meas-tech.net/5/457/2012/ 


A. J. M. Piters et al.: The CINDI campaign: design, execution and early results 


477 


3648-element ID CCD was used to record the spectra. Its 
temperature was kept constant at 40.0±0.1°C. Integration 
time was kept constant throughout the day at around 100 ms. 
Spectra recorded at a fixed elevation angle for a 5 min in- 
terval were averaged and analysed. The telescope unit was 
placed on the roof of one of the temporary cabins. A 
single planoconvex lens (focal length = 40 mm and diame- 
ter = 25 mm) was used to focus the received sunlight onto the 
fibre. Window and lens are both made of quartz without op- 
tical coatings. A movable mirror (UV reflecting) in the tele- 
scope unit was controlled by a PC to record spectra sequen- 
tially at six different elevation angles. Power consumption 
was roughly estimated to be less than 200 VA for the whole 
system, including a PC. 

NASA Pandora 2, 3, and 9 

Pandora’s head sensor consists of a quartz entrance win- 
dow, one or two 9-position filter-wheels (opaque, UV-optical 
band-pass filters, polarizers, neutral density filters), and a 
quartz lens focusing on a single-strand fibre. The head sen- 
sor is mounted on a tracker capable to point at any position 
in the sky up to a zenith angle of 1 10°. The 10 m long fibre 
with NA of 0.12 acts as a depolarizer and transmits the light 
to the Avantes mini-spectrometer. Pandora has non-cooled 
Hamamatsu detectors (CCD or CMOS). The spectrometer 
is temperature stabilized to 10 ± 3 °C. The complete system 
weighs about 12 kg. The maximum power consumption is 
200 W. Pandora 9 is a sun-only instrument. It took direct-sun 
data from sunrise to sunset using different optical band-pass 
filters in the filter-wheel. A direct-sun measurement consists 
of multiple spectra at typically 10 ms integration time for a 
total of 40 s. Pandora 2 and 3 took 1 direct-sun measurement 
per 5 elevation scans. 

NIWA MAX-DOAS 

The NIWA MAX-DOAS instrument consists of an 
ISA HR320 (a 0.32 m Czerny-Turner monochromator) 
fitted with additional internal baffles. A Hamamatsu 7420 
detector with S7032-1007 cooled — 20 °C back-thinned 
CCD with 1044 x 128 pixels is used with FFT-CCD in line 
binning mode. The dark signal is le - /pixels -1 and the 
readout noise 8e“s _1 . The telescope optics consists of 
a fibre optic with multiple 100 pm diameter fibres, which 
removes polarization, and an F/4 lens. The fibre optic con- 
nects vertically to an elevation-azimuth tracker. The power 
usage of the NIWA MAX-DOAS was 120 W. The instrument 
was installed in a metal box of somewhat less than lm 3 . 
The instrument box had a thermoelectric cooler and was 
additionally cooled with room air from the cabin during the 
latter part of the campaign when ambient temperatures were 
high. 


Universities of Leicester and Leeds MAX-DOAS 

The MAX-DOAS operated jointly by the Universities of Le- 
icester and Leeds during the CINDI campaign briefly com- 
prised of a 25 mm diameter telescope, built at the University 
of Leeds, with a 7.5 cm focal length mounted on a stepper 
motor for elevation angle control. The telescope and step- 
per motor were housed within a heated protective case with 
the telescope connected via a 10 m fibre optic bundle to a 
temperature stabilized PI Acton Insight Spectrometer inside 
the building. The spectrometer, with a grating of 1200 lines 
per mm, was housed in a case stabilized to ±0.1K. The 
CCD detector used (PI Acton model: Pixis 1024) is back il- 
luminated and contains 1024 rows and 1024 columns with 
a pixel size of 13 x 13 pm 2 giving a total imaging area of 
13.3 x 13.3 mm 2 . The camera is kept in a permanent vacuum 
and maintained at a constant temperature of —42 °C using a 
thermo-electric Peltier cooler. Spectra were collected from 
averages over the central 500 rows of the CCD. 

University of Toronto Ground-Based Spectrometer 
(UT-GBS) 

The University of Toronto Ground-Based Spectrometer (UT- 
GBS) collects sunlight with a 40 mm-diameter fused silica 
lens (f-number 2.5). The light is focused onto aim liq- 
uid light guide, with a core diameter of 3 mm, which de- 
polarizes the light and brings it to a three-lens f-number 
matcher. This focuses light onto a 0.1 mm slit and through 
a mechanical shutter into an Instruments SA/Jobin Yvon 
Triax-180 crossed Czerny-Turner triple grating spectrome- 
ter (Fraser et al., 2009). The back-illuminated CCD detector 
with 2048 x 5 12 pixels is thermoelectrically cooled to 20 1 K. 
The instrument operated in two viewing modes - zenith sky 
(ZS) and multiple elevation angles (MD) - and two differ- 
ent wavelength ranges (UV and VIS). For the MD mode a 
custom-built sun-tracker was placed above the instrument. 
During twilight, only zenith measurements were taken. The 
UT-GBS was housed on the ground in an air-conditioned alu- 
minium box (69 x 126 x 67 cm 3 ), weighing 100 kg total. To- 
tal power consumption was 2000 W of power at 120 V. 

Washington State University MFDOAS instrument 

MFDOAS instrument (WSU) measures direct-sun (DS) ir- 
radiance and scattered sky (SS) radiance using an imag- 
ing Czerny-Turner spectrometer (Acton SpectraPro 2356, 
focal length 300 mm) and a back illuminated, UV-coated 
charge coupled detector (Princeton Instruments PIXIS- 
2KBUV, 512 x 2048 pixels 2 (6.9 x 27.6 mm 2 ). SS pho- 
tons are collected by an off-axis parabolic mirror telescope. 
The DS field of view is defined by a 2.5 cm biconvex 
lens and baffles inside the tube connected to the integrat- 
ing sphere. Before entering the spectrometer, light passes 
through two filter wheels that contain UV transmitting filters 
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(Hoya U340) for removing visible light and depolarizers 
(Halbo Optics WDQ25).The observation scheme included 
vertical scanning combined with direct-sun measurements 
after 23 June. The instrument positioner (Kipp and Zo- 
nen 2APG) is capable of 360° azimuthal and 90° zenith ro- 
tations. The CCD was cooled to —70 °C and the instrument 
enclosure was temperature regulated at 28± 1 °C. The instru- 
ment was located outside at the RSS. The instrument dimen- 
sion is 0.8 x 1 x 1.8m 3 , with a total weight ~100kg. The 
power consumption is 1.58 kW. 

IUP-Heidelberg MAX-DOAS 1 and 2 

The MAX-DOAS instrument 1 of IUP Heidelberg consists 
of a separate telescope and spectrometer unit, which are con- 
nected by a quartz fibre bundle. It was installed at the wind 
profiler site, where the telescope unit was set-up outdoors 
and the temperature controlled spectrometer unit indoors un- 
der stable conditions. The incoming scattered sunlight is 
collected via two quartz prisms, whose rotation axes are or- 
thogonal, allowing for the pointing to any direction in the 
sky. The light is focused by a spherical mirror and dis- 
tributed to the three miniature Ocean Optics HR2000 spec- 
trometers with a symmetrical crossed Czerny-Turner optical 
bench with a f/# of 4 and a linear 2048 pixel CCD array. 
The MAX-DOAS instrument 2 has a Czerny-Turner spectro- 
graph (Acton 300). Light was detected using a 2-dimensional 
back-illuminated CCD array with 2048 x 256 pixels. View- 
ing direction sequences with an elevation angle of 2, 4, 8, 15, 
30 and 90° were repeatedly performed, with an integration 
time per viewing direction of one minute. 

A2 SAOZ instruments 

CNRS SAOZ 

The standard SAOZ (Pommereau and Goutail, 1988) as in- 
stalled at the Wind Profiler site is made of a commercial 
Jobin-Yvon CP200 flat field spectrometer equipped with a 
holographic grating of 360 grooves mm -1 associated to a 
non-cooled 1024 NMOS diode array detector from Hama- 
matsu and a 50 pm entrance slit. The SAOZ is placed in a 
sealed tight and dehydrated box on the top of which a quartz 
window is mounted to view at zenith sky. The system is com- 
pletely automated. A GPS is used for time and solar zenith 
angle (SZA) accurate calculation. The exposure time is ad- 
justed automatically between 0.1 and 60s in order to opti- 
mize the signal and the spectra are co-added in memory dur- 
ing a 60 s duty cycle. During the campaign, measurements 
were performed every 2 min up to a SZA of 94°. The tech- 
nical specifications are the following: 70 cm diameter, 40 cm 
height, 28 kg weight and a power consumption of 500 W at 
220 V. 


CNRS mini-SAOZ 

The zenith-sky mini-SAOZ, the new miniaturized version of 
SAOZ and the first prototype, also participated in the cam- 
paign. It is made of a commercial Czerny-Turner spectrome- 
ter from Avantes with a 2048 x 14 non-cooled CCD detector. 
Measurements were performed every 2 min. The technical 
specifications are the following: 20 x 15 cm 2 , 10 cm height, 
2 kg weight and 300 W power at 220 V. 

CNRS mobile SAOZ 

The mobile SAOZ is the balloon version of SAOZ powered 
by batteries (in order to allow measurements onboard a car) 
and in which the conical optical head used usually for direct- 
sun observations was replaced by a quartz window and a field 
of view limiter of 9° for zenith-sky observations. The scien- 
tific specifications are identical to the standard SAOZ. The 
housing is different and the technical specifications are the 
following: 50 x 45 cm 2 , 90 cm height, 20 kg weight. 

A3 In-situ instruments 

EMPA in-situ 

EMPA measured NO and NO 2 in-situ with a modified com- 
mercial chemiluminescence analyser TEI 42C TL (Thermo 
Fischer Scientific Inc.). Measurements were taken continu- 
ously from 12 June to 3 July 2009 at 10 min time resolution. 
Air was sampled at an altitude of about 3 m above ground and 
drawn to the instrument installed in the basement of the tower 
through a Teflon tubing of about 10 m length. The instrument 
measures NO by detecting the chemiluminescence generated 
in the reaction of NO with excess ozone, while NO 2 is first 
reduced to NO by means of a converter (Kley and McFar- 
land, 1980). For the selective measurement ofN02, the stan- 
dard molybdenum converter was replaced by a photolytic 
(blue light) converter (Droplet Measurement Technologies, 
Boulder, Co, USA). Automated daily calibrations included 
zero and span cycles using a reference gas (20 ppm NO in 
N 2 , Messer-Griesheim, Germany) diluted with purified (NO- 
free) ambient air, and measurements of the conversion effi- 
ciency by gas phase titration. The conversion efficiency was 
of the order of 45 % during the whole campaign. The total 
uncertainty (1 <r) of the 10-min samples c were estimated to 
V 7 0. 13 2 -(- (0.035c) 2 ppbv and yj 0.29 2 + (0.078c) 2 ppbv for 
NO and NO 2 respectively. 

RIVM and LML in-situ 

During the CINDI campaign, RIVM operated three chemi- 
luminescence nitrogen oxides analysers. All three were 
Model 200E from Teledyne Instruments. One analyser was 
part of the operational Dutch National Air Quality Monitor- 
ing Network (LML), and operated in its standard configu- 
ration with a molybdenum converter to convert NO 2 to NO. 
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The other two were equipped with photolytic converters, also 
from Teledyne Instruments. The monitor with molybdenum 
converter was located in the base of the tower at Cabauw, and 
had its air inlet at 3 m above the ground. One of the monitors 
with photolytic converter was using the same air inlet for the 
first half of the campaign. After 10 July 2009 this monitor 
was moved to the 100 m level of the tower. The other mon- 
itor with photolytic converter was located at the 200 m level 
of the tower for the entire campaign. 

IUP-Bremen in-situ 

IUP Bremen operated a commercial nitrogen oxide analyser 
(Eco Physics CLD 88p) with a photolytic NO 2 converter 
from the same company (PLC 860). During the campaign 
the same calibration source as for the EMPA instrument has 
been used. The conversion efficiency was of the order of 
60 % throughout the campaign. The Bremen instrument was 
located in the basement of the tower from 12 June 2009 to 
21 July 2009 with the inlet (Teflon tubing) at 3 m altitude 
outside the building. 

A4 Lidars 

RIVM Caeli lidar 

The lidar Caeli (CESAR water vapour, AErosol and cloud 
Lidar, Apituley et al., 2009) is a high-performance, multi- 
wavelength Raman lidar, capable of providing round-the- 
clock measurements. The instrument is part of the Eu- 
ropean Aerosol Research Lidar Network (EARLINET), 
and provides profiles of volume backscatter and extinc- 
tion coefficients of aerosol particles, the depolarization ra- 
tio, and water-vapour-to-dry-air mixing ratio. A high-power 
Nd:YAG laser transmits pulses at 355, 532, and 1064 nm. 
Because a large telescope is essentially blind for lidar signals 
from close to the instrument, a second, small telescope is 
used to cover the near range, in particular for measurements 
in the planetary boundary layer. 

RIVM mobile NO 2 lidar 

The RIVM mobile lidar instrument uses the DIAL technique 
(Gimmestad, 2005) to measure NO 2 concentrations. The in- 
strument uses a dye laser to send laser pulses, alternating be- 
tween 413.463 nm and 414.112nm, into the atmosphere. A 
280 mm telescope, narrow band interference filter and photo- 
multiplier tube comprise the detector. The laser beam can be 
pointed in any direction. The measurement range is typically 
from 300 m to 2.5 km, sometimes up to 4 km. During the 
CINDI campaign, measurements at elevation angles of 0.5°, 
1°, 2°, 4°, 8°, 15°, 30°, and 90° were combined into a single 
concentration profile. Such a profile takes 5 minutes or more 
to measure, runs from 5 m to up to 4 km, and has a vertical 
resolution ranging from several meters at the bottom to 2 km 
at the top of the profile. The entire system is housed in a 


fully self-supporting mobile laboratory, 8 m long, 2.5 m wide 
and 2.3 m high, mounted on a vehicle. During the CINDI 
campaign, it was located on the south-western tip of the Re- 
mote Sensing Site. The instrument is extensively described 
by Volten et al. (2009). 

RIVM tropospheric ozone lidar 

The Tropospheric Ozone Lidar, Trop03, was built and op- 
erated at RIVM in the Netherlands in the 1990’s (Sunesson 
et al., 1994) for the purpose of routine monitoring of tro- 
pospheric ozone. Trop03 is a differential absorption lidar 
(DIAL) based on two frequency-quadrupled Nd:YAG lasers 
(FHG), emitting 266 nm laser light. Stimulated Raman scat- 
tering in pressurized gas cells is used to create the DIAL 
wavelengths /, (m at 289 nm and /. 0 ir at 299 nm. The backscat- 
tered light is received by a 60 cm diameter telescope from 
which ozone profiles are retrieved between approximately 
2 km altitude and 15 km altitude. A second receiver with a 
diameter of 5 cm is used to measure ozone profiles between 
about 500 m and 2 km. A full tropospheric ozone profile is 
measured by averaging data during about 20 to 30 min. Con- 
ditions with broken cloud cover of up to 5 octa can be used 
to provide ozone profiles by rejection of cloud contaminated 
data. The spatial resolution of the retrieved ozone profiles 
varies between 100 m at low altitudes and 1 km towards the 
top of the troposphere. 

A5 Sondes 

KNMI N0 2 sondes 

Experimental NO 2 sondes have flown on 6 days: 18, 23, 
24, 25, 30 June, and 1 July, all launched at Cabauw around 
10:30 UTC. The sonde uses the chemiluminescent reaction 
of NO 2 in an aqueous luminol solution, which is optimised 
to be specific to NO 2 (Sluis et al., 2010). The sonde is at- 
tached to a small meteorological balloon. It has a vertical 
resolution of 5 m and a measurement range between 1 and 
lOOppbv. The instrument weighs 0.7 kg. 

KNMI O 3 sondes 

ECC ozone sondes were launched on the same days as the 
NO 2 sondes at 12:30 UTC. The ECC ozone sensor (Komhyr, 
1969; Komhyr and Harris, 1971) is an electrochemical cell 
consisting of two half cells, made of Teflon, which serve 
as cathode and anode chamber, respectively. Both half cells 
contain a platinum mesh serving as electrodes. They are im- 
mersed in KI solution of different concentrations. The two 
chambers are linked together by an ion bridge in order to 
provide an ion pathway and to prevent mixing of the cathode 
and anode electrolytes. The ECC does not require an external 
electrical potential. The ECC gets its driving electromotive 
force from the difference in the concentration of the KI solu- 
tion in the cathode and anode chamber. The electrical current 


www.atnios-meas-tech.net/5/457/2012/ 


Atmos. Meas. Tech., 5, 457-485, 2012 


480 


A. J. M. Piters et al.: The CINDI campaign: design, execution and early results 


is directly related to the uptake rate of ozone. The sonde is 
flown in a polystyrene protective box (source: Harris et al., 
1998). 

Appendix B 

Meteorological parameters 

Daily meteorological values during CINDI are listed in Ta- 
ble Bl. For each day the following parameters are listed: av- 
erage cloud cover, total number of sun hours, average wind 
direction and wind speed, maximum temperature, average 
pressure, hours of precipitation, average relative humid- 
ity, and average NO 2 surface volume mixing ratio between 
04:00 and 19:00 UTC (day time). All parameters were mea- 
sured in Cabauw. The cloud-cover is the average between 
06:00 and 18:00 UTC, derived from TSI data (Long et al., 
2006). The average NO 2 was determined from the in-situ in- 
struments from RIVM, IUP Bremen, and EMPA. Days with 
at least 10 sunny hours are defined as “category A” days, days 
with at least 4 sunny hours are defined as “category B” days 
(indicated in column “cat”). A sunny hour is defined as an 
hour with less than 50% cloud cover and with more than 
0.5 h of sunshine, as derived from the hourly meteorological 
values measured in Cabauw. 

Appendix C 

Satellite overpasses 

Table Cl lists the satellite overpasses during CINDI for 
GOME-2 aboard EUMETSAT METOP, for SCIAMACHY 
aboard ESA ENVISAT and for OMI aboard NASA EOS- 
Aura. Category A and B days (Sect. 3.1) are indicated 
in the second column. At some days GOME-2 and OMI 
had a second overpass. These are not indicated in the ta- 
ble. The retrieved tropospheric NO 2 column and its un- 
certainty are given for the pixel closest to Cabauw. These 
values are obtained via the Tropospheric Emission Moni- 
toring Internet Service (http://www.temis.nl). No overpass 
time means: no measurement within 50 km, no NO 2 value 
means: retrieval not trusted because of clouds (cloud radi- 
ance fraction >50%, i.e. cloud cover >~20%) or the so- 
called “row-anomaly” (for OMI, indicated by an asterisk). 
Processor versions are: TM4N02A version 2. 10 for GOME- 
2, TM4N02A version 2.0 for SCIAMACHY, DOMINO ver- 
sion 2.0 for OMI. 
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Table Bl. Average meteorological values during CINDI, see description in Appendix B. 


Date 

Cat 

Cc 

(%) 

Sun 

(h) 

Wind 

dir 

Wind 

speed 

(ms -1 ) 

Temp 

(C) 

Press 

(hPa) 

Precip 

(h) 

Rel 

hum 

(%) 

no 2 

(ppbv) 

8 June 2009 

B 

61 

9.6 

182 

3.3 

18.4 

1006.9 

0.9 

73 

8.7 

9 June 2009 


82 

4.6 

212 

5.3 

19.6 

1005.2 

4.9 

82 

8.2 

10 June 2009 


99 

1.0 

223 

3.0 

17.6 

1011.1 

6.4 

85 

9.8 

1 1 June 2009 


83 

2.8 

302 

4.1 

17.0 

1013.7 

8.4 

85 

2.5 

12 June 2009 

A 

13 

14.4 

298 

3.7 

17.6 

1021.0 

0.0 

73 

2.9 

13 June 2009 

A 

20 

12.7 

179 

1.8 

21.6 

1021.1 

0.0 

69 

3.9 

14 June 2009 


93 

4.0 

250 

3.0 

21.0 

1017.7 

4.6 

84 

4.8 

15 June 2009 


93 

3.3 

7 

1.7 

18.9 

1016.6 

1.3 

80 

6.5 

16 June 2009 


73 

5.2 

5 

2.6 

18.4 

1021.4 

1.1 

82 

7.2 

17 June 2009 

B 

14 

10.2 

174 

3.5 

23.1 

1019.8 

0.0 

71 

7.5 

18 June 2009 

B 

41 

7.8 

269 

4.0 

20.4 

1016.7 

0.4 

69 

4.1 

19 June 2009 

B 

53 

9.5 

271 

5.1 

18.1 

1018.5 

0.0 

70 

2.7 

20 June 2009 


98 

6.7 

284 

4.1 

17.5 

1021.4 

0.0 

75 

2.1 

21 June 2009 

B 

93 

8.4 

338 

3.4 

17.7 

1022.8 

0.5 

77 

1.7 

22 June 2009 

B 

63 

10.4 

6 

2.5 

19.0 

1027.6 

0.0 

74 

3.4 

23 June 2009 

A 

22 

14.8 

21 

4.3 

21.3 

1027.7 

0.0 

70 

4.7 

24 June 2009 

A 

38 

14.3 

23 

5.1 

22.8 

1020.9 

0.0 

72 

5.0 

25 June 2009 

B 

41 

9.1 

64 

3.7 

26.3 

1013.7 

0.0 

64 

5.7 

26 June 2009 


53 

9.6 

64 

2.9 

24.8 

1011.7 

0.0 

77 

4.2 

27 June 2009 


98 

4.4 

332 

1.8 

25.4 

1014.4 

0.0 

84 

6.3 

28 June 2009 


93 

2.2 

318 

1.8 

24.2 

1016.0 

0.0 

91 

3.6 

29 June 2009 

B 

71 

8.8 

329 

2.5 

26.1 

1018.8 

0.0 

85 

3.5 

30 June 2009 

A 

22 

13.4 

360 

3.0 

25.2 

1021.4 

0.0 

82 

5.1 

1 July 2009 

B 

38 

9.9 

19 

3.2 

26.2 

1022.2 

0.0 

80 

4.1 

2 July 2009 

A 

37 

12.8 

48 

3.0 

30.2 

1017.9 

0.0 

71 

5.4 

3 July 2009 

B 

52 

10.2 

240 

4.3 

28.5 

1012.9 

0.0 

74 

8.4 

4 July 2009 

A 

8 

14.2 

306 

2.5 

25.2 

1013.2 

0.0 

72 

3.5 

5 July 2009 


75 

8.4 

200 

2.2 

27.7 

1010.3 

0.4 

72 

4.2 

6 July 2009 


71 

8.9 

231 

4.0 

23.1 

1007.3 

0.5 

77 

4.5 

7 July 2009 


95 

2.8 

214 

5.5 

18.2 

1003.9 

4.9 

86 

5.7 

8 July 2009 


78 

6.8 

282 

5.8 

19.1 

1008.3 

0.5 

80 

2.3 

9 July 2009 


84 

4.5 

266 

4.0 

18.5 

1013.0 

0.0 

75 

3.8 

10 July 2009 


92 

2.1 

280 

5.2 

17.8 

1013.3 

6.4 

82 

3.0 

1 1 July 2009 


81 

5.7 

256 

3.2 

20.3 

1013.8 

0.1 

78 

5.2 

12 July 2009 


73 

5.6 

216 

4.4 

21.4 

1009.0 

6.6 

86 

4.5 

13 July 2009 


59 

11.3 

196 

3.4 

23.1 

1010.8 

0.2 

73 

6.1 

14 July 2009 

B 

57 

8.4 

237 

2.0 

24.4 

1010.7 

0.0 

73 

7.5 

15 July 2009 

B 

38 

10.9 

222 

4.8 

24.2 

1015.9 

0.0 

71 

5.1 

16 July 2009 

A 

27 

13.9 

170 

3.1 

25.7 

1018.6 

0.0 

65 

6.3 

17 July 2009 


69 

6.9 

202 

5.0 

22.7 

1007.2 

3.0 

72 

8.0 

18 July 2009 


83 

4.5 

237 

6.0 

19.6 

1008.4 

4.1 

82 

3.3 

19 July 2009 


74 

7.3 

232 

5.3 

20.1 

1011.9 

1.5 

77 

2.4 

20 July 2009 

B 

54 

8.8 

242 

5.1 

21.0 

1014.7 

1.3 

75 

4.5 

21 July 2009 


85 

6.5 

148 

3.2 

26.8 

1009.7 

1.4 

80 

10.5 

22 July 2009 


80 

4.8 

248 

3.1 

24.2 

1004.6 

5.1 

84 

9.4 

23 July 2009 


89 

2.9 

226 

2.8 

20.5 

1005.0 

3.4 

86 

9.0 

24 July 2009 


70 

7.5 

230 

5.0 

20.3 

1011.8 

1.2 

80 

3.9 

25 July 2009 

A 

28 

13.5 

273 

4.3 

20.6 

1020.5 

0.4 

79 

3.4 

26 July 2009 


74 

8.0 

195 

2.8 

23.4 

1018.9 

0.0 

74 

5.4 
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Table Cl. Satellite overpasses during CINDI, see description in Appendix C. 


Date 

Cat 

GOME-2 

SCIAMACHY 

OMI 

Time 

[UTC] 

N0 2 (err) 
[10 15 cm -2 ] 

Time 

[UTC] 

N0 2 (err) 
[10 15 cm“ 2 ] 

Time N0 2 (err) 

[UTC] [10 15 cm -2 ] 

8 June 2009 

B 

09:47 




11:55 


9 June 2009 


09:26 


10:39 


12:38 

8(2) 

10 June 2009 


10:45 


10:08 


11:43 


1 1 June 2009 


10:24 




12:26 


12 June 2009 

A 

10:04 

9(4) 

10:44 

11(7) 

11:31 

4.4 (1.9) 

13 June 2009 

A 

9:43 

7(5) 

10:14 

12(9) 

12:13 

6(4) 

14 June 2009 


9:22 




11:19 


15 June 2009 


10:41 




13:39 


16 June 2009 


10:21 




12:44 

2(2) 

17 June 2009 

B 

10:00 

16(10) 



11:49 

14(4) 

18 June 2009 

B 

09:39 




12:31 

16(11) 

19 June 2009 

B 

09:19 


10:25 

3.0(1. 5) 

11:37 


20 June 2009 


10:38 




12:19 


21 June 2009 

B 

10:17 




11:25 


22 June 2009 

B 

09:57 




12:07 

3.4(1. 8) 

23 June 2009 

A 

09:36 

5(5) 

10:00 


12:49 

* 

24 June 2009 

A 

09:15 




11:55 

5(3) 

25 June 2009 

B 

10:35 

12(6) 

10:37 

8(5) 

12:38 

6(4) 

26 June 2009 


10:14 

12(5) 

10:06 

17(7) 

11:43 


27 June 2009 


09:53 




12:25 

9(6) 

28 June 2009 


09:32 


10:42 


11:31 


29 June 2009 

B 

09:12 


10:10 


12:13 

13(11) 

30 June 2009 

A 

10:31 

5(3) 



11:19 

8(2) 

1 July 2009 

B 

10:11 

15(9) 

10:49 

8(4) 

13:39 

5(3) 

2 July 2009 

A 

09:50 

5(3) 

10:17 

6(4) 

12:44 

7(4) 

3 July 2009 

B 

09:29 




11:49 


4 July 2009 

A 



10:53 

8(5) 

12:31 

5(4) 

5 July 2009 


10:28 


10:23 

6(4) 

11:37 


6 July 2009 


10:07 




12:19 

6(4) 

7 July 2009 


09:46 




11:25 


8 July 2009 


09:26 


10:28 


12:07 


9 July 2009 


10:45 




12:49 

* 

10 July 2009 


10:24 




11:55 


1 1 July 2009 


10:04 


10:35 


12:37 

4(3) 

12 July 2009 


09:43 


10:02 


11:43 


13 July 2009 


09:22 




14:04 

4(3) 

14 July 2009 

B 

10:41 


10:40 

13(8) 

11:31 


15 July 2009 

B 

10:21 


10:08 


13:51 

5(2) 

16 July 2009 

A 

10:00 

14(9) 



11:18 

10(8) 

17 July 2009 


09:39 


10:44 


13:39 


18 July 2009 


09:19 


10:14 


12:43 

1.3 (1.8) 

19 July 2009 


10:38 




11:49 


20 July 2009 

B 

10:17 

4.8 (1.6) 

10:51 

6(3) 

12:31 

6(3) 

21 July 2009 


09:57 




11:37 


22 July 2009 


09:36 




12:19 

8(6) 

23 July 2009 


09:15 




11:24 


24 July 2009 


10:34 


10:25 


12:07 


25 July 2009 

A 

10:14 




12:49 

* 

26 July 2009 


09:53 

8(5) 



11:55 



* indicates that there is an OMI overpass, but the measurement is affected by the so-called “row-anomaly”. 
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